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NATIONAL ADVISOR^ CCKklTZte FOR AERONAUTICS 
TECHNICAL NOTE NO . 1370 

CORREIATION OF EXEERIMElNTAL AND CALCULATED EFFECTS OF 
PRODUCT OF INERTIA ON LATERAL STABILITY 


By Marlon 0. McKinney, Jr. and Hubert M. Drake •. . 


SUMMARY 


A correlation of experimental and calculated effects of the 
product of inertia on the stability of idle lateral oscillations of 
airplanes has been obtained from flight teste of a free -flying air- 
plane model having a 42® sweptback wing and from calculations of the 
stability of the model. In order to provide a comprehensive correla- 
tion, the static directional stability of the model was varied by 
changing the vertical -tail area and tail length and the inclination 
of idle longitudinal principal axis of inertia relative to the wind 
axis was varied by changing the wing incidence. 

The calculated lateral -stability boundaries were found to be 
in good agreement with measured lateral stability when idle product- 
of -inertia terms were included in the calculations. Neglecting the 
product of inertia, however, led to wide discrepancies between 
calculated and measured stability. These results emphasized the 
necessity for considering the product of inertia in lateral -stability- 
analyses. 

The general flying characteristics were influenced primarily 
by the static directional stability and to a lesser degree by the 
stability of the lateral oscillations. A certain minimum amount 
of static directional stability was required to give good controlla- 
bility even when less directional stability provided good damping 
of the lateral oscillations. 


INTRODUCTION 

Until recently the effects of the product of inertia . usually 
have been neglected in lateral -stability analyses because the 
lateral -stability studies of reference 1 indicated that these 
effects were relatively unimportant for conventional airplanes. 
Calculations in reference 2, however, show that the product of 
inertia may have a pronounced effect on the lateral stability of 
scone high-speed airplanes because of high wing loadings, large 
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differences ‘between yatfirig arid roiling moments of inertia, high 
operational altitudes, ana sweepback . The sweepback may cause 
high effective dihedral and high angles of attack and, consequently, 
large angles between the principal longitudinal axis of Inertia 
and the wind axis. 

In order to obtain an experimental cheolc of lateral -stability 
calculations including the product -of -inertia terms, a systematic 
series of flight tests have been made in the Langley free -flight 
tunnel. A model having a 42° sweptback wing was used to determine 
experimentally the effects of the product of inertia on the lateral 
stability of the free -flying model for correlation with the cal- 
culated stability characteristics of the model. In order to 
provide, a comprehensive check of the calculations , the directional 
stability of the model was varied, by changing the vertical-tail 
size and tail le ngth a nd the quantitative effects of the product 
of inertia were varied by changing the wing incidence and thereby 
changing the inclination of the principal axes of inertia relative to 
the wind axes. 


SYMBOLS 


The forces and moments are referred to the stability axes, which 
are defined as an orthogonal system of axes intersecting at the 
airplane center of gravity in which the Z-axis is in. the plane of 
symmetry and perpendicular^ to the relat ive wind, the X-axis is in 
the. plane of symmetry and perpendicular* to the Z-axis, and the 
Y-axis is perpendicular to the plane of symmetry. A diagram of 
these axes showing the positive direction of forces and moments 
is presented in figure 1. 

The symbols and coefficients are defined as follows: 

W . weight of model, pounds 

m mass of model, slugs 

S wing area, square feet 

S^. vertical -tail area, square feet 

b wing span, feet 

c wing chord, feet 

tail length (distance from center of gravity to 
rudder hinge line), feet 


l 
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Sr 

% 

% 

%Z 

Y 

q. 

p 

r 

P 

P 

a 

r 

n 

C L 

Cy 

C, 


n 

°Yi 


0 


height of center of pressure of vertical tail above 
fuselage axis, feet 

angle of incidence of ving, degrees 

radius of gyration of model about principal longitudinal 
axis, feet 

radius of gyration of model about principal normal axis, 
feet 

product -of -Inertia factor, feet 2 ((*Z 2 - % 2 ) cos T) sin t}) 

angle of attack of principal longitudinal axis of airplane; 
positive vhen forward end of major principal axis Is 
above X-axis 


(H> 


airspeed, ■ feet- per second 
dynamic pressure, pounds per square foot 
rolling angular velocity, radians per second 
yawing angular velocity, radians per second 
mass density of air, slugs per cubic foot 
angle of sideslip, degrees except -where otherwise noted 
angle of attack of fuselage axis, degrees 
angle of climb, degrees 
airplane rolative -density factor 


lift coefficient 


(if) 

lateral -for ce coefficient 


vW 

\ _ qS 


) 


rolling -moment coefficient 


( 


Bolling moment 


q£>b 


) 


( yawjyiA VKtoQQ £ 


yawing -moment coefficient ^ 


5 


rate of change of lateral -force coefficient with angle of 
sideslip in radians 
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C, 


n ( 


C Yr 
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rate of change of rolling -moment coefficient with angle 
of Bideslip in radians (except where specified to he 

in degrees) 

rate of change of yawing -moment coefficient with angle of 
sideslip in radians (except where specified to he in 

degrees) 

rate of change of lateral -force coefficient, with rolling - 

/acO 

angular -velocity factor in radians I 

\ 2V t 

rate of change of rolling -moment coefficient with rolling - 
angular-velocity factor in radians 




rate of change of yawing -moment coefficient, with rolling- 
angular -velocity factor in radians 



rate of change of rolling -moment coefficient with yawing- 
angular-velocity factor in radians 



rate of change of yawing-moment coefficient with yawing - 

_ /ac n \ 

angular -velocity factor in radians j — if f 

YW 


APPARATUS 


The investigation was conducted in the Langley free-flight tunnel, 
which is equipped for testing free -flying dynamic airplane models* A 
description of the tunnel and its operation is given in references 3 
and 4 . Free -oscillation tests were made to determine the damping-in- 
yaw derivative C„ "by the method described in reference 5 * Steady- 

rotation tests were made to determine the damping -in -roll deriva- 
tive by the method described in reference o . 

A sketch of the model used in the investigation is shown ae 
figure 2. Ihis model is described in reference 7 except for a few 
changes which wore necessary for the present investigation. Wedge- 
shape blocks were provided in order that the wing might be mounted 
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on the fuselage with 0°, 10°, and -10° incidence. The nodel was 
equipped with a set of interchangeable vertical tails, ranging in 
size from 2.6 to 10.4 percent of the wing area, which could be 
mounted at any of the three tail-length positions indicated in 
figure 2. The various model configurations tested are identified 
in table I. 


TESTS AND CALCULATIONS 


Force tests of. the model were made to determine the values of 
static -lateral -stability derivatives with various vertical-tail 
arrangements and wing incidences. Flight tests of the model were 
made at a lift coefficient of 0.6 with angles of incidence of the 
wing of 0°, 10°, and -10° and with vertical -tail arrangements 
given in table I that provide the values of the directional - 
stability parameter CL. and the effective dihedral parameter -C 7 

3 ■ 

indicated in figure 3* 


£ 


The model was flown at each test condition to determine the 
stability of the Dutch roll oscillation, which was recorded by 
photographing the uncontrolled motions of the model after a dis- 
turbance caused by abruptly deflecting the ailerons to roll the 
model from approximately 30 ° bank to level flight and then centering 
the controls. In some of the cases - the flying characteristics of 
the model were too poor to permit recording the uncontrolled motions 
of the model and it was necessary to resort to the pilot. r s observa- 
tions to determine whether the model was stable. The pilot's 
ingress ion of the general flying characteristics of the model were 
also noted in order to provide data concerning the ease of flying 
the model, both for level flight and for performance of the mild 
maneuvers possible in the tunnel. The ratings of the general flying 
characteristics of the model in controlled flight generally indicate 
vhother stability and controllability are adequate and properly 
proportioned. The vertical -tail area and tail length were selected 
in an attempt to bracket the oscillatory-stability boundaries in 
order to provide the desired correlation of test and theory. 

All the flight tests were made with the rudder coupled to the 
ailerons for lateral control, and the ratio of rudder deflection to 
aileron defleqtion was adjusted for each test condition to minimize 
the adverse yawing due to rolling and aileron deflection. For the 
lower values of directional stability, however, the adverse yawing 
could not be entirely eliminated. 
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Calculations wore mad© "by the method presented in reference 2 
to determine the "boundary of neut./al stability of the lateral oscilla- 
tion of the model. Two sets of calculations were made, one taking 
into consideration and the other neglecting the product of inertia. 

The boundaries for the model with a 10° angle of incidence of the 
wing were computed by assuming changes of the tail size, and the 
boundaries for the model with 0° and -10° angles of incidence of 
the wing were computed by assuming changes in the tail length in 
order to make the calculations consistent with the test procedure. 


(tall off) were obtained from 
for configurations B and C 


The aerodynamic and mass characteristics of the model used in 
the calculations are presented in table II. The mass characteristics 
of the model were obtained by measurements. The inclination of the 
principal longitudinal axis of inertia to the fuselage axis was 
found to be less than 0.3° for all v-Ihq incidences and was assumed 
to be zero in the calculations. The trim airspeed, flight-path angle, 
and angle of attack were determined from free -flight tests. The 
values of Cy (tail off) and C nn 
p p 

force tests . The values of Cy„ , 

^(tail) 

were also obtained from force tests. (See table II.) The tail -off 
values of C^ and C n _^ were determined from steady- rotation and 

free -os dilation tests, respectively; and the tail -off values 
of C* r and were estimated from the charts of reference 8 and 

unpublished wind-tunnel data. The tail contributions to the stability 
derivatives were estimated from the equations given in the footnote 
of table II, which are similar to those given in reference 9, by 
utilizing the experimentally determined values of Oy contributed 

3 

by the vertical tails and the arbitrarily chosen values of tail 
length in the case of configurations B and C and by utilizing the 
given tail length and arbitrarily chosen values of Cy in 

P (tail) 

the case of configuration A. 


RESULTS AND DISCUSSION 
Stability 


The results of the stability investigation are presented in 
figures 3 and 4. Figure 3 shows the calculated oscillatory -stability 
boundaries of the model and indicates whether the uncontrolled 
lateral oscillations of the model were stable, neutral, or unstable 
in the configurations tested. Figure 4 shows typical time histories of 
the uncontrolled rolling motion of the model following a rolling 
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disturbance . Only the rolling motions sore presented for simplicity. 
The period and damping of the rolling motions,, however, are identical, 
with those of the yawing and sideslipping. Inasmuch as spiral 
stability is considered relatively unimportant end the calculations 
indicate no effect of the products of inertia, the spiral stability 
of the model was not investigated. 

■When the product of inertia terms were included in the calcula- 
tions of the stability boundaries, the boundaries were in good 
agreement with the flight test results as shown in figures 3 (a) 
and 3(b). Figure 3(c) shows that the model was stable, as is 
indicated by the calculated stability boundary which included 
product of inertia. However, it was impossible to check this, 
boundary closely because the model could not be flown at negative 
values of C n . A simple illustration of the pronounced effect of 
3 ■ 

the product of inertia may be obtained by a comparison of test 
points for configurations and By in which the model con- 

figurations were identical except for the change in wing incidence 
and the consequent change in the angle of attack of the principal 
longitudinal axiB of inertia. In configuration. A 3 (fy. = 10°) 
the model was unstable whereas in configuration By - 0°) the 

model was quite stable. 

..The test results are at variance with the stability boundaries 
which were calculated by neglecting the product of inertia. This 
difference is particularly evident in figures 3(h) and 3(c). Some 
additional calculations were made in order to insure that poor esti- 
mates of the stability derivatives which were not measured could not 
account for the discrepancy between the test results and the calcu- 
lations in which the product of inertia wan neglected. The results 
of the calculations are presented in figure 5 which shows the effects 
on the stability boundary of introducing a large value of the 
stability derivative Cy^, doubling the tail -off value of , 

or using one -half the estimated value of Cn (tail off). None of 

■these changes in the stability derivatives could account for more 
than approximately one -fourth of the difference between the calcula- 
tions neglecting the product of inertia terms and the point of 
neutral stability as determined from the tests. The product of 
inertia terms . should therefore be included in lateral-stability • 
calculations . 

The stability or instability of the uncontrolled motions of the 
models, as shown in figure 4, was quite definite except for test 
configuration B 3 which appears, to be about neutrally stable. No 
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records of the uncontrolled lateral motions could "be made for the 
condition of -10° wing incidence because the poor flying charac- 
teristics of the model with so littlo directional stability made 
such records very difficult to obtain. The pilot's impression 
was, however, that the lateral oscillations of the model for these 
conditions were heavily damped. 

The data of figures :3 and 4 are in good agreement with refer- 
ence 2 regarding the pronoijinced effects of the product of inertia 
on lateral stability. At the present time there are not sufficient 
data to determine when the effects of the . product of ihox-tia .will 
be important^ hence; the product of inertia should be considered 
in all lateral -stability analyses. 


General Flight Behavior 

The flight -test- results. presented in figure 6 show that the 
general flight behavior of- the model was influenced primarily by 
the static directional stability, which has a pronounced effect 
on the response of. the model to the controls, and to a lesser degree 
by the dynamic lateral stability. 

At constant C^, increasing the oscillatory stability by 

decreasing the wing incidence improved the general flight behavior 
of the model because an increase in the stability increases the 
tendency of the model to fly itself and thereby mahes the piloting 
of the model easier and the flight- smoother. The most apparent 
effect of the oscillatory stability on the general flight behavior 
was found by comparing the flying characteristics of configurations A 
and B-j_ 

Increasing the directional stability improved the general 
behavior of the model by improving the response to controls as 
was explained in reference 10 . This effect was especially noticeable 
with the present model because of the high effective dihedral of the 
model. When adverse yawing occurred this high effective dihedral 
produced rolling moments which, opposed the aileron rolling moments 
and thus reduced the effectiveness of the ailerons for controlling 
the model. The data presented in figure 6 show that the model was 
easier to fly in configurations A2 and than in configurations B2 
and B3, in spite of the fact that the lateral oscillations were 
unstable for the A configurations and stable or neutrally stable 
for the B configurations. The directional stability of the model 
was too low for satisfactory response to the controls in ccnfigura- . 
tions B2 and B3. For configuration C3 the response of the model 
to the controls was so poor that it was virtually unflyable although 
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the lateral oscillations were apparently very stable, as was Indicated 
by the calculations . The model could not be flown with any less 
directional stability than it had in configuration Cg. 

Xt was therefore concluded that general flight behavior was 
influenced primarily by static directional stability and to a lesser 
degree by the dynamic lateral stability and that a certain min i mu m 
amount of static directional stability was required to give good 
response to the controls even when less directional stability 
provided good damping of the lateral oscillations. 


.CONCLUSIONS 


The following conclusions were drawn from an investigation in • 
the Langley free -flight tunnel to obtain a correlation between the 
calculated and -the measured effects of the product of inertia on 
lateral stability characteristics; 

1. The calculated- lateral -s tability boundaries were in good 
agreement with measured lateral stability when the product, of 
inertia terms were included in the calculations * Neglecting the 
product of inertia, 'however, led to wide discrepancies between 
calculated and measured stability. These results emphasized the 
necessity for considering the product of inertia in lateral- 
stability analyses. 

2. The general flying characteristics were influenced primarily 
by the static directional stability and to a lesser degree by the 
stability of the lateral oscillations . A certain minimum amo unt 

of static directional stability was required to give good controlla- 
bility even when less directional stability provided good damping of 
the lateral oscillations. 

Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., May 21, 1947 
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•VERTICAL-1A.IL AERAlT02 /t m , JTS FOR THE VARIOUS 
MODEL COimGURATEOMS 


Model 

configuration 

Wing 

incidence, 

. 

(deg) 

Vertical- 
tail area. 

st/s 

Tail 

length, 

1/b 

A 1 

10 

0.104 

0.67 

Ag 

10 

• 

o 

CD 

.67 

A 3 

10 

.052 

■67 

B 1 

0 

.052 

.67 

B 2 

0 

.05a 

.41 

B 3 

0 

.052 

•15 

C 1 

-10 

.026 

.67 

C 2 

-10 

.026 

.41 

°3 

-10 

.02 6 

•15 


R ATIO NAL ADVISORY 
COl'f'HTIEE FOR AEROHAUECCS 
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tabu n 

CHARACTERISTICS 07 TEB MDXZLOSED ID BI OAICDIATIOHS 



Configuration A 

Configuration B 

Configuration C 

V 

ID .0 

O 

-10.0 

W 

5.10 

b.6b 

5.10 

v/s 

2.0b 

1.855 

2.0b 

b 

3.6E 

3-82 

3.82 

P 

o .00238 

0.00238 

0.00238 

V 

53-5 

51.2 

53-5 

p. 

7-00 

-.34 

7.00 

% 

0.V75 

0.510 

0 -b 75 

V 

1-315 

1.290 

1-315 

*xz 

O.Obb7 

*0.191 

-0-b38 

1/b 

0.®T 

Variable 

Variable 

*/b 

0 . i 3 hi x FcZl 

V 3 (tall) 

0.072 

0.036 

°L 

0.6 

0.6 

0.6 

a 

-b.0 

6.0 

16.0 

7 

-9.0 

- 9.0 

-9.0 

% 

^ * %««« 

-0.007b + Cr 

^P(tall) 

-0.007b + C, 

~P(tail) 

*% 

° + °V(tail) 

0 + C- 

0 + 0- 

P(tall) 


-0.230 + C,_ 

T’(tall) 

- 0.230 + C, 

V(tall) 

-0.230 + Ci 

S(tall) 

*% 

■°-° 33S+ % ( teil) 

-0.033a + 0 

^(tall) 

-0.0332 + c 

jP(tall) 

a ° l r 

0.125 + C, 

0.125 + C, 

lr (tall) 

0.125 + Cr 

r (tail) 


-0.0090 + c 

-0.0090 + c_ 

^■(tail) 

-0.0090 + C„ 

^■(tail) 

016 (tall) 

Variable 

-0.1870 

-0.078b 


"tell contribution* arc determined from the relieving equations: 

l NATIONAL ADVISORY 

S(tail) “ "V %(tall) COHKITTEE FOt AERONAUTICS 


C, m £(*■ - ala a) G~ 

V(tall) V * > ^S(tail) 

°V(tail) " °* r (tall) ’ ^ * Bl “ O^Oftall) 

^(tall) °^®(tall) 

Varied eyetematioally as Independent variable to provide the desired range of Ojj. for the determination 
of tho b tab ill t^r ■bcomdaarf.eB * p 
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Fig. 1 


Y 



Figure /. - 7he stability system of axes. Arrows indicate 
posi/ire directions of moments, forces , and control- 
surface deflections. This system of axes is defined as 
an orthogonal system haying fhe/r origin at the 
center of grar/ty ond in which the Z-axts is /n the 
plane of symmetry and perpendicular to the relative 
wind , the X - axis is in the plane of symmetry and 
pierpendicular to the Z - axis > and the Y~ axis is 
perpendicular fo the p/one of symmetry 
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Fig. 3 



£ ffective-dihedral parameter, <\ p ) per deg 

Figure 3— Correlation of oscillatory-stability ratings from flight tests with calculated oscillatory- 

stability boundaries. 
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Figured- Flight records of the uncontrolled rolling motions of the model 
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Fig. 5 


Oscillatory-stability ratings ^ O Ota 6/e 
from flight tests ^ Neutral 



NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 

Figure S. — Effect of varying some of the stability 
derivatives an the oscillatory -stability boundaries 
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(sta6k) Calculated osciHatorv-stabthtv boundaries 

(Urstabk) ' ' including produc ts of inertia 
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